US  Army  Corps 
of  Engineers 


Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to 
the  originator. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other 
authorized  documents. 


The  contents  of  this  report  are  not  to  be  used  for 
advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of 
such  commercial  products. 


Preparation  of  this  report  was  sponsored  by  the  Department  of  the  Army 
In-House  Laboratory  Independent  Research  Program.  Neither  the  Depart¬ 
ment  of  the  Army,  nor  the  US  Army  Corps  of  Engineers,  nor  their  employees, 
nor  any  of  their  contractors,  subcontractors,  or  their  employees  make  any 
warranty,  expressed  or  implied,  or  assume  any  legal  liability  or  responsibility 
for  the  accuracy,  completeness,  or  usefulness  of  any  information,  apparatus, 
product,  or  process  disclosed,  or  represent  that  its  use  would  not  infringe 

privately  owned  rights. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704^188 

PubtK  repoain^  burden  for  This  cofieaiorr  of  information  is  estimated  to  aver*^  l  hour  oer  resoonse.  including  the  time  for  reviewing  instruaions.  searching  existing  data  sources, 
gathcnng  ar>d  maintainirtg  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  asoect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Heaoguaners  Services.  Directorate  for  information  Operations  and  Reports.  I2  IS  iefferson 
Oavis  Highway.  Suite  1204.  Arlington,  VA  22202-4)02.  and  to  the  Office  of  Management  and  Oudget.  Paperwork  Reduction  Project  (0704-0180).  Washington.  DC  2050). 

1.  AGENCY  USE  ONLY  (Leivt  bisnk)  2.  REPORT  OATE  _  3.  REPORT  TYPE  AND  OATES  COVERED 

September  1992  Report  1  of  a  Series 

4.  TITLE,  AND  SUBTITLE 

Application  of  the  Green-Naghdi  Theory  of  Fluid 

Sheets  to  Shallow-Water  Wave  Problems;  Report  1, 

Model  Development 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Zeki  Demirbilek,  William  C.  Webster 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

USAE  Waterways  Experiment  Station,  Coastal  Engineering 
Research  Center,  3909  Halls  Ferry  Road,  Vicksburg,  MS 
39180-6199 

Ship  Research,  Inc.,  55  Highland  Blvd. ,  Kensington,  CA 
94707 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

Technical  Report 
CERC-92-11 

9.  SPONSORING/ MONITORING  AGENCY  NAME(S)  AND  AOORESS<ES) 

US  Army  Corps  of  Engineers,  Washington  DC  20314-100 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

Available  from  National  Technical  Information  Service, 
Springfield,  VA  22161. 

5285  Port  Royal  Road, 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  worcb) 

This  report  presents  the  mathematical  formulation 
a  numerical  model  that  simulates  wave  transformation  ir 
model  is  intended  both  for  military  and  civil  worl 
propagation  of  time -dependent  and  nonlinear  waves  wher 
either  be  inapplicable  or  use  of  simple  analytic  or  r 
infeasible.  The  theory  detailed  in  this  report  introduces 
wave  model  for  shallow  to  moderate  water  depths  where  the 
The  Creen-Naghdi  Level  II  theory,  hereafter  referred  to  s 
has  been  significantly  modified  in  this  research  and  a  pov 
numerical  model,  called  CNWave,  is  developed  for  water  w 
The  theory  and  ensuing  model  incorporate  some  o 
mathematical  features  of  the  water  wave  equations, 
approximating  of  the  governing  Euler's  field  equations  a 
boundary  conditions  necessary  for  capturing  the  bulk  phys: 

for  the  development  of 

1  shallow  waters .  The 
cs  projects  involving 
e  existing  models  may 
lumerical  solutions  is 
a  new- generation  water 
seabed  varies  rapidly, 
imply  as  the  CN  theory, 
rerful,  general-purpose 
ave  problems . 

:  the  most  important 
These  include  non- 
nd  imposing  the  proper 

Lcal  characteristics  of 

/  CnnM  rmpdl _ 

14.  SUBJECT  TERMS 

See  reverse 

»  _ _ _ _ _ _ 

IS.  NUMBER  OF  PAGES 

48 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

NSN  7540-01 -280-5500  Standard  Form  298  (Rev  2-89) 


Pf«cnbed  by  ANSI  Std  i39*i8 


290-102 


13.  (Concluded). 


wave  trains  in  the  shallow- water  regime.  The  GN  approach,  which  is  fundamentally 
different  from  the  perturbation  method  based  on  developments  in  classical  wave 
theory  begun  by  Stokes  and  Boussinesq  in  the  last  century,  can  do  this  only 
because  it  does  not  introduce  any  simplifications  of  the  velocity  variation  in 
the  vertical  direction  across  the  fluid  layers  or  sheets. 

In  contrast  to  the  Stokes  and  Boussinesq  theories,  the  equations  of  motion 
in  the  GN  theory  are  derived  by  enforcing  exact  kinematic  and  d3rnamic  boundary 
conditions  on  the  free  surface  and  on  the  bottom,  and  by  enforcing  conservation 
of  mass  and  of  the  0th  and  1st  moments  of  momentum  in  the  vertical  direction. 
These  conditions  yield  11  coupled  partial  differential  equations,  which  can  be 
reduced  to  3  complicated  governing  equations  by  elimination  of  many  of  the 
variables.  In  summary,  the  GN  theory  is  different  from  the  perturbation  approach 
in  that  the  free  surface  and  bottom  boundary  conditions  are  met  exactly,  whereas 
the  field  equation  is  implicitly  approximated.  The  result  is  a  theory  that  can 
predict  the  shape  and  behavior  of  waves  up  to  almost  breaking  conditions.  The 
GN  theory  breaks  down  when  the  particle  velocity  at  the  crest  equals  the  wave 
speed,  the  criterion  for  breaking  in  the  exact  theory. 

By  developing  the  unrestricted  Green-Naghdi  theory  of  fluid  sheets ,  this 
research  presents  a  new  wave  theory  consisting  of  a  coupled,  nonlinear  set  of 
partial  differential  equations  and  integrates  these  in  time  and  space  to  simulate 
either  regular  or  irregular  real  waves.  The  theory  and  model  have  been  shown  to 
reproduce  with  engineering  accuracy  the  evolution  of  a  wave  of  permanent  form, 
from  small  amplitudes  up  to  almost  breaking  conditions.  The  theory  presented  is 
for  a  nonlinear  numerical  wave  tank  in  which  the  seabed  topography  profile  can 
be  arbitrary  and  very  irregular,  and  up  to  20  wave  gages  can  be  positioned  at 
will  inside  the  computational  domain  to  obtain  snapshots  and  profiles  of  wave 
records , 

The  types  of  projects  to  which  the  theory  can  be  applied  are  many,  and 
include  problems  of  both  military  and  civil  interest.  The  theory  is  purposely 
made  to  be  versatile  to  permit  decision  makers,  designers,  and  analysts  to  assess 
the  various  aspects  of  waves  and  wave -structure  interaction  problems  arising  in 
Army  applications.  One  can  evaluate,  for  instance,  the  effect  of  submerged 
obstacles  during  military  landings  on  the  train  of  waves  approaching  a  beach  or 
landing  zone,  or  the  reflection  of  waves  and  forces  on  sea  walls  or  spillway 
hydraulic  gates,  and  the  time  history  of  bottom-mounted  pressure  gage 
measurements  for  estimation  of  surface  wave  conditions  in  coastal  design 
projects.  The  theory  is  particularly  suited  for  the  violent  collision  of  waves 
with  natural  and  man-made  structures,  and  their  impact  on  preventive  and 
defensive  hydraulic  structures. 
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APPLICATION  OF  THE  r.RHEN-NAGHDl  THEORY  OF  FLUID  SHEETS  TO 
SHALLOW-WATER  WAVE  PROBLEMS 


PART  I:  THEORETICAL  BASIS 
Introduction 

This  is  the  first  report  in  a  series  prepared  for  the  Technical  Director,  U.S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  for  the  project  entitled  Time-Dependent  Nonlinear  Free-Surface 
Wave  Modeling  conducted  under  the  Director’s  auspices  as  part  of  the  In-House  Laboratory 
Independent  Research  Program  (Project  A91D  per  Office  of  Technical  Planning  and  Programs 
guidelines).  This  report  is  a  joint  effort  between  researchers  at  the  WES  Coastal  Engineering 
Research  Center  (CERC)  and  the  University  of  California,  Berkeley,  under  Purchase  Order 
DACA39-91-M-5515.  It  is  intended  to  provide  a  discussion  and  critique  of  various  approaches  for 
simplifying  complex  hydrodynamic  boundary  value  problems,  a  derivation  of  the  general 
Green-Naghdi  (GN)  theory  of  fluid  sheets,  theoretical  documentation  of  the  equations  of  motion  of 
two-dimensional  (2-D)  shilow-water  waves  for  GN  levels  1  and  2,  and  a  concise  description  of  the 
numerical  methods  used  to  integrate  these  governing  equations. 

The  intent  of  the  authors  was  to  assemble  a  document  that  a  fluid  mechanician  can  use  to 
understand  this  relatively  new  method  and  its  application  to  complex  shallow-water  wave  problems. 
The  entire  theoretical  formulation  presented  in  this  report  is  new,  greatly  improving  upon  several 
previously  published  papers  or  Ph.D.  dissertations  completed  at  the  University  of  California, 

Berkeley.  The  derivation  of  the  theory  is  in  the  form  of  a  tutorial  in  which  all  of  the  intermediate 
^>teps  are  included,  since  no  textbook  or  article  is  available  in  this  level  of  detail. 

During  the  two  decades  since  its  introduction  (Green,  Laws,  and  Naghdi  1974),  the  theory  of 
fluid  sheets  has  been  applied  to  a  variety  of  fluid  flow  problems.  These  include  studies  of  waves  in 
shallow  and  deep  water  (Green  and  Naghdi  1986  and  1987;  Shields  and  Webster  1988),  the  flow 
beneath  planing  boats  (Naghdi  and  Rubin  1981),  the  waves  created  by  a  moving  pressure  disturbance 
(Ertekin  1984;  Ertekin  et  al.  1984  and  1986)  and  wave  reflection  by  obstacles  (Marshall  and  Naghdi 
1990),  to  name  a  few.  In  particular,  the  development  of  fluid  sheet  theory  in  a  Eulerian  frame 
(Green  and  Naghdi  1984)  made  this  theory  much  easier  to  apply  to  fluid  flow  problems.  The  reader 
is  referred  to  the  pair  of  papers  by  Green  and  Naghdi  (1986  and  1987)  for  a  definitive  exposition  of 
the  theory,  and  especially  for  its  application  to  water  wave  problems. 

Overview  of  Theory 

Alternative  approaches,  with  a  variety  of  approximations  and  assumptions,  exist  for 
calculating  wave  motion  in  coastal  waters.  TTie  classical  equations  of  motion  for  fluid  flow  in  three 
dimensions  are  a  continuum  model  which  embodies  many  assumptions.  For  ordinary  fluids,  such  as 
water,  the  Navier-Stokes  equations  are  universally  recognized  as  a  good  model  for  the  resulting  flows. 
However,  these  equations  are  not  "exact"  equations  but  are  an  idealization  similar  in  spirit  to  the 
idealization  of  space  by  Euclidean  geometry.  Even  for  simple  free-surface  problems,  these  equations 
and  their  simpler  inviscid  counterparts,  the  Euler  equations,  are  difficult  to  solve.  One  popular 
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approach  has  been  to  systematically  simplify  the  three-dimensional  equations  and  their  boundary 
conditions  through  a  formal  perturbation  analysis  until  the  resulting  system  can  be  solved.  The 
theories  of  water  waves  developed  by  Stokes,  Boussinesq,  and  others  follow  this  type  of  development. 


Green-Naghdi  fluid  sheet  theory  offers  an  alternative  in  the  form  of  a  new  model,  that  of  a  2- 
D  continuum  of  unsteady  3-D  flows.  Although  the  examples  which  will  be  cited  here  involve  inviscid 
fluids,  the  development  of  GN  theory  is  not  at  all  limited  to  such  fluids.  The  discussion  below  is 
aimed  at  exploring  the  difference  between  these  two  very  different  paths  to  simplification  of  the 
analysis  of  fluid  flow  problems.  In  either  case,  it  is  anticipated  that  the  solutions  obtained  are 
approximate  ones,  since  there  really  is  no  substitute  for  solving  the  three-dimensional  (3-D)  equations 
exactly.  Both  approaches  are  called  approximations,  although  it  is  clear  that  the  meaning  is  not  the 
same  for  each. 

Before  introducing  details  on  the  nature  of  GN  theory,  it  is  useful  to  first  discuss  the  notion  of 
approximation  in  general.  An  approximation  approach  for  analyzing  a  given  problem  is  u  ally 
chosen  based  on  its  ability  to  predict  the  phenomena  that  one  is  interested  in  and  on  its  ease  of  use. 
The  selection  of  an  approximation  scheme  can  be  viewed  as  a  type  of  non-zero-sum  game  where  one 
attempts  to  i.'ake  assumptions  which  will  have  a  greater  impact  on  the  simplification  of  the  analysis 
than  on  the  accuracy  of  the  prediction  of  the  phenomena  of  interest. 

Two  observations  from  this  discussion  are  significant.  First,  the  choice  of  the  approximation 
scheme  depends  on  the  specific  answers  that  one  is  looking  for  (i.e.,  the  choice  depends  on  the 
context  of  the  problem  rather  then  its  generic  type).  Second,  the  means  of  analysis  change  in  time; 
that  is,  computations  which  20  years  ago  would  have  required  the  world’s  largest  computers  can  now 
be  accomplished  faster  and  for  a  minimal  cost  on  a  personal  computer.  It  is  proper  to  think  in  terms 
of  approximation  schemes  "appropriate  for  the  current  time."  Since  the  evolution  of  a  new  computer 
generation  appears  to  take  only  a  few  years,  it  seems  natural  that  we  will  see  a  corresponding 
evolution  in  approximation  schemes  which  will  take  advantage  of  these  new  resources.  It  is  a  thesis 
of  this  report  that  GN  theory  and,  in  particular,  higher-level  GN  theory  is  appropriate  for  our  time. 
Two  developments  lead  to  this  conclusion;  the  emergence  of  low-cost,  high-speed  computation,  and 
the  emergence  of  sophisticated  symbolic  manipulation  software  which  allows  one  to  accurately 
perform  calculus  and  algebraic  manipulations  on  rather  large  systems  of  equations. 

Approximation  schemes  can  be  separated  into  different  categories.  Perturbation  methods,  both 
ordinary  and  singular,  introduce  some  mathematical  approximation  to  reduce  the  complexity  of  the 
model  to  the  point  where  it  can  be  solved.  One  advantage  of  these  methods  is  that  one  obtains 
governing  equations  for  the  flow  and  from  these,  both  specific  solutions  can  be  obtained  and  general¬ 
izations  of  the  behavior  of  the  flow  can  be  made.  On  the  other  end  of  the  spectrum,  the  original 
problem  can  be  solved  by  purely  numerical  techniques.  Finite  difference,  finite  element,  and  panel 
methods  are  such  schemes.  These  methods  are  comparable  to  physical  experiments  in  that  each 
computation  yields  another  result  corresponding  to  a  single  realization  of  the  flow.  Generalization 
about  the  behavior  of  the  flow  requires  induction  from  many  of  these  specific  solutions.  GN  fluid 
sheet  theory  lies  in  the  middle  of  this  spectrum.  It  achieves  simplification  by  reducing  the 
dimensionality  from  three  dimensions  to  two.  This  theory  yields  governing  equations  for  the  flow, 
which  are  solved  numerically  in  a  more  efficient  manner  than  those  from  the  three-dimensional 
model. 
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In  perturbation  analyses,  reference  scales  appropriate  for  the  particular  problem  at  hand  are 
introduced.  These  scales  are  used  to  nondimensionalize  the  variables  and  to  identify  a  nondi- 
mensional  perturbation  parameter  (or  parameters)  which  can  be  considered  small  (or  large).  For  time 
invariant  problems,  the  flow  is  decomposed  into  a  sequence  of  flows  of  presumably  decreasing  impor¬ 
tance,  each  of  which  is  a  correction  to  the  sum  of  the  previously  computed  flows.  The  assumed 
sequence  is  inserted  into  the  field  equations  and  boundary  conditions  and  the  perturbation  parameters 
are  used  to  segregate  these  into  a  corresponding  sequence  of  perturbation  problems.  Typically,  each 
of  these  problems  is  linear  in  the  unknowns  at  its  level,  although  it  may  involve  higher-order  terms  of 
quantities  determined  already  in  previous  (lower-order)  solutions. 

An  implicit  assumption  is  made  that  this  sequence  is  convergent,  but  this  is  almost  never 
proven.  In  some  flow  problems,  such  as  two-dimensional  water  waves  in  both  shallow  and  deep 
water,  there  is  ample  evidence  of  the  convergence  (Schwartz  1974).  In  problems  such  as  the  flow 
about  thin  airfoils,  the  lack  of  convergence  is  well  known.  These  methods,  which  date  back  to  Stokes 
and  other  early  researchers,  are  often  called  "rational  methods"  because  the  assumptions  are  clear  and 
testable,  and  the  details  can  be  embodied  into  a  mathematical  process  through  which  theoretically  one 
can  obtain  solutions  to  whatever  level  of  accuracy  one  chooses,  if  the  perturbation  sequence 
converges.  A  particular  advantage  of  the  perturbation  approach  is  that,  since  perturbation  parameters 
are  used  to  size  quantities,  the  ingredients  of  this  parameter  give  one  insight  into  the  types  of 
problems  for  which  the  approximation  is  appropriate.  However,  the  perturbation  approach  does  not 
yield  quantitative  measures  of  the  accuracy  to  be  expected  for  a  particular  problem.  This  information 
can  only  be  obtained  from  an  analysis  of  higher  order  problems  or  from  comparison  with 
experiments. 

Unsteady  fluid  problems  are  rather  different.  It  is  usually  not  feasible  to  consider  the  flow  as 
a  sum  of  linear  perturbation  problems  that  one  can  solve  sequentially  until  sufficient  accuracy  has 
been  obtained,  unless  some  additional  limiting  assumption  such  as  periodic  motion  is  introduced. 
Generally,  one  must  solve  a  single  set  of  governing  equations  in  time.  Perturbation  methods  have 
also  been  used,  for  instance,  by  Wu  (1981)  for  the  formulation  of  approximations  appropriate  for 
time-domain  wave  problems.  For  these  problems,  the  introduction  of  scales  permits  grouping  of 
terms  of  like  size  for  a  particular  problem.  One  can  obtain  a  variety  of  different  sets  of  governing 
equations  depending  on  the  order  of  terms  retained.  That  is,  one  can  obtain  a  sequence  of  sets  of 
governing  equations  for  a  given  problem,  each  of  which  contains  all  of  the  terms  of  the  previous  sets 
plus  those  due  to  the  retention  of  the  next  order  of  smaller  terms.  Presumably,  this  sequence  of 
increasing  complexity  will  produce  solutions  of  increasing  accuracy.  Generally,  all  of  these  sets  of 
governing  equations  will  be  nonlinear  with  the  exception  perhaps  of  the  first.  However,  since  we  are 
throwing  away  parts  of  the  exact  problem,  we  can  expect  that  some  quantities,  such  as  mass  and 
momentum,  may  not  exactly  be  conserved  (although,  if  a  consistent  analysis  has  been  performed,  the 
eifors  should  be  of  a  size  comparable  to  the  first  neglected  order). 

The  limitations  resulting  from  such  an  analysis  can  be  subtle.  Many  perturbation  schemes 
consider  the  fluid  velocity  to  be  a  small  perturbation  to  a  reference  velocity.  The  typic-tl  result  is  that 
any  order  of  the  perturbation  theory  is  not  Galilean  invariant,  since  the  terms  which  are  needed  to 
make  it  so  are  spread  amongst  several  orders.  One  such  example  is  the  Korteweg  de  Vries  equations 
for  nonlinear,  shallow-water  flow. 

In  problems  where  viscosity  is  not  important  (or  can  be  ignored)  and  the  flow  is  initially 
quiescent,  the  field  equation  becomes  Laplace’s  equation,  which  is  linear.  For  these  flows  much  of 
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the  focus  of  perturbation  analysis  is  therefore  on  the  boundary  conditions.  Both  the  kinematic 
boundary  condition  on  material  surfaces,  as  well  as  the  dynamic  boundary  condition  on  the  free 
surface,  are  nonlinear.  Expansions  of  these  nonlinear  conditions  and  grouping  of  terms  by  orders  of 
the  perturbation  parameter(s)  lead  to  a  basis  for  selecting  or  discarding  terms.  For  these  flows,  then, 
the  perturbation  method  solves  a  problem  where  the  exact  field  equation  is  satisfied,  but  where  the 
boundary  conditions  are  satisfied  only  approximately.  In  some  sense,  the  uniqueness  of  these  flow 
situations  stems  solely  from  the  imposed  conditions  at  the  boundaries  and  it  is  worrisome  to  concen¬ 
trate  the  approximations  there. 

Green-Naghdi  theory  is,  however,  quite  the  opposite  in  nature.  The  boundary  conditions  are 
met  exactly,  but  the  field  equations  are  approximated.  In  this  approach  the  dependence  of  the 
kinematic  structure  of  the  solutions  along  one  coordinate  direction  is  prescribed.  This  direction  is  the 
vertical  (or  x^)  direction  for  the  theories  discussed  in  this  paper.  In  many  problems,  such  as  shallow 
water  problems,  the  depth  of  the  fluid  in  the  x^  direction  may  be  quite  small;  in  others  the  fluid 
domain  can  be  infinite  in  this  direction.  In  either  case  the  resulting  theories  are  called  fluid  sheet 
theories.  The  assumed  variation  of  the  fluid  velocity  across  x’  will  be  expressed  here  as  a  finite  sum 
of  products.  The  first  term  in  each  product  is  a  coefficient  that  depends  on  the  remaining  two 
horizontal  coordinates  (x'  and  x^  and  time,  and  the  second  term  is  a  function  of  x’  alone. 

Approach 

The  governing  equations  for  the  Green-Naghdi  theories  presented  here  are  composed  of;  an 
exact  statement  of  the  conservation  of  mass,  an  approximate  statement  of  the  conservation  of 
momentum,  and  exact  statements  for  the  various  boundary  conditions.  As  a  result  of  the  formulation, 
x’  no  longer  appears  in  the  governing  equations  and  all  quantities  therein  are  functions  of  x',  x%  and 
time.  No  scales  are  introduced  and  no  terms  are  thrown  out.  Development  of  GN  theory  takes  place 
in  two  steps:  postulation  of  a  set  of  governing  equations,  and  verification  that  these  equations  satisfy 
certain  physical  requirements.  The  process  presented  below  is  one  which  is  unique  to  this  project, 
and  although  it  closely  follows  research  efforts  in  the  subject  area  at  the  University  of  California, 
Berkeley,  it  is  substantially  different  from  that  used  by  Green  and  Naghdi. 

The  first  step  is  a  procedure  for  identifying  a  candidate  fluid  sheet  model.  In  fact,  there  is  no 
preferred  method  for  this  identification  and  at  this  point  the  model  could  just  as  well  have  been 
induced  from  the  results  of  model  tests.  The  same  variational  approach  of  Kantorovich  and  Krylov 
(1958)  used  previously  by  Shields  and  Webster  (1988)  is  used  to  derive  the  candidate  model  from 
standard  three-dimensional  equations.  This  approach  is  a  variation  of  the  method  of  weighted 
residuals,  and  is  therefore  similar  in  nature  to  the  procedure  used  in  the  development  of  finite 
elements.  In  this  procedure,  the  dimensionality  of  the  system  of  partial  differential  equations  is 
reduced,  rather  than  the  system  being  replaced  by  the  system  of  algebraic  equations,  as  it  would  be  in 
a  Galerkin  procedure. 


PART  II:  FORMULATION  WITH  GENERAL  WEIGHT  FUNCTIONS 


Governing  Equations 

Let  X  =  x'  (i  =1,2,3)  be  a  system  of  fixed  Cartesian  coordinates  in  Euclidian  space  with  base 
vectors  Cj,  where  63  is  oriented  vertically  upward.  For  convenience,  x^  is  denoted  by  f  in  the 
subsequent  development  because  this  dimension  plays  a  much  different  role  than  the  other  two 
dimensions.  In  the  following  standard  Cartesian  tensor  notation  is  used,  with  the  summation 
convention  implied  for  repeated  indices.  In  many  instances,  however,  the  summation  will  be  stated 
explicitly  for  clarity.  Latin  indices  are  used  for  quantities  having  three  spatial  components  and  take 
on  values  of  1,  2,  3;  Greek  indices  take  on  the  values  of  1  and  2  only.  A  comma  in  the  subscript 
denotes  differentiation  by  the  following  variable  or  that  corresponding  to  the  subsequent  index. 

The  fluid  velocity  vector  at  a  point  x  and  time  t  is  given  by  v  =  v  (x,  t)  =  V;  Cj.  The  fluid  is 
assumed  to  be  bounded  by  two  smooth  and  non-intersecting  material  surfaces.  The  material  surfaces 
are  given  by  f  =  a(x‘,  x^,  t)  and  f  =  /3(x‘,  x^,  t),  /3  >  a,  respectively.  Since  a  and  /3  are  material 
surfaces,  on  these  surfaces  the  kinematic  ("no  leak”)  boundary  condition  is 

Da  r  ,, 

—  =  [V3  -  a,  - 

f  =  h  -  p.  -  ',P,llc., 

This  discussion  will  be  concerned  only  with  an  incompressible  ideal  fluid,  and  in  this  case  the 
stress  vector  t  =  -p  e;  and  the  mass  density  of  the  fluid  p  are  assumed  to  be  constant.  The  body  force 
is  given  by  -  pge3.  (In  a  more  general  theory,  t  would  have  the  form  tj  =  t  y  Uj,  where  T|j  is  the 
stress  tensor,  and  the  body  force  would  be  a  general  vector  pf.)  Unit  outward  normal  vectors  on  the 
top  and  the  bottom  surfaces  are  denoted  by  A  and  n,  respectively.  The  three-dimensional  (Euler) 
equations  resulting  from  the  conservation  laws  of  mass  and  momentum  are 


pv^  +  (pv,v)^  =  -p^e,  -  pgej 

The  fundamental  kinematic  assumption  that  the  velocity  field  can  be  approximated  is 
introduced  as 


K 

v(Jt‘,ac*,C,0  =  5^  W,(x^x^^)  A.,(0 

fi«0 


where 

W,(x‘,x^r)  =  Wl  (.x\x\t) 

X„((')  is  a  "shape"  function  that  depends  upon  f  only.  The  coefficients  W„  are  unknown  time- 
dependent  vectors  to  be  determined  as  a  part  of  the  solution.  The  W„  correspond  with  the  "directors" 
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in  the  original  work  (Green  and  Naghdi  1986  and  1987).  For  each  choice  of  K  a  complete,  closed  set 
of  equations  is  developed  that  is  independent  from  those  for  a  different  value  of  K.  Thus,  the 
kinematic  models  form  a  hierarchy  depending  on  K  and  increasing  in  complexity  with  K.  Since  this 
hierarchy  is  different  from  a  perturbation  expansion,  the  terminology  adopted  is  that  suggested  later 
(Shields  and  Webster  1988)  to  describe  the  complexity  of  the  theory.  A  particular  member  of  this 
hierarchy  is  referred  to  as  the  "K"*  level  approximation." 


The  kinematic  boundary  conditions  in  Equation  1  may  be  rewritten  using  Equation  3  as 


Y:  w'AM)  =  a,  -  E 

ll»0  AsO 

E  -  p,  •  E 


(4) 


A«0 


The  continuity  equation.  Equation  2a,  likewise  becomes 

r  jr 


E  K,  1.  *  E  ■  0 


(5) 


In  addition  to  the  kinematic  boundary  conditions,  4K  scalar  equations  are  needed.  K  scalar 
equations  derived  from  conservation  of  mass  are  chosen  and  K  vector  equations  derived  from 
conservation  of  momentum  (corresponding  to  3K  scalar  equations  for  three-dimensional  problems)  to 
provide  closure  for  the  fluid  sheet  theory.  It  is  convenient  at  this  point  to  restrict  the  weighting 
functions  to  those  that  possess  the  following  property 

^m.i  =  E  ^ 

r*0 


where  the  a,  are  constants.  The  function  set  is  therefore  a  finite  closed  set  under  differentiation. 
Inserting  Equation  6  into  Equation  5,  the  continuity  equation  can  be  expressed  as 

E  I  "  E  <  ^.(C)  =  0  (7) 

(  »i»0  J 

or,  since  the  terms  in  braces  are  not  a  function  of  f 

"  E  <  °r|  =  0  (8) 


for  r  =  1,...,  K  . 
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Equation  8  is  therefore  an  exact  statement  of  conservation  of  mass  for  the  flow  given  by  the 
kinematic  approximation,  Equation  3.  Note  that  if  the  derivative  of  the  r*  weighting  function  with 
respect  to  f  is  not  expressible  in  terms  of  the  previous  orders  of  the  weighting  function  set,  more  than 
K  conditions  will  result  from  this  procedure  and  the  Krylov-Kantorovich  method  described  below 
could  be  used  to  determine  approximate  equations  for  the  conservation  of  mass. 

There  are  many  function  sets  that  satisfy  Equation  6  (for  instance,  exponential  functions 
where  m  =  n).  Polynomial  functions  also  satisfy  Equation  6,  like  those  previously  used  for  shallow- 
water  fluid  sheet  theory.  Similarly  the  sets  {sinh(af),  cosh(af)}  and  {sin(ai'),  cos(af)  },  a  =  ao,  ai,  ... 
,  a„  also  satisfy  Equation  6. 

If  Equation  3  is  substituted  into  momentum  Equation  2b,  and  the  resulting  equation  is 
required  to  be  satisfied  everywhere  in  the  fluid  domain  (as  for  continuity),  many  more  equations  than 
the  desired  K  vector  equations  would  be  obtained.  This  difficulty  is  due  to  the  presence  of  the 
quadratic  terms  in  Equation  2b.  A  weak  formulation  due  to  Kantorovich  is  employed  where  the 
"shape  functions"  X„  are  used  as  weighting  functions  to  develop  K  approximate  equations  which 
express  the  conservation  of  momentum  in  some  integral  sense.  Multiplying  Equation  2b  by  each  Xn(f) 
and  integrating  through  the  vertical  direction  results  in 

//  [(P  V),  (p  V,  V),]  A.(0^c  =  [  -P/,  -  pge,]  X„(0^/C  w 


for  n  =  1,...,K. 

Using  the  product  rule  of  differentiation,  and  noticing  that  \  is  not  a  function  of  t,  and  that  y 
ranges  from  1  to  2,  Equation  9  can  be  expressed  as 

/J  [(  P  v)^  =  (pvjv)  ,.]  X,  dC 

"  I!  a  ^  I! 

which  further  reduces  to 

It  is  noted  that  X„  is  only  a  function  of  f  and  this  means  that  X„,.,  =  0.  With  this  relation  and 
recognizing  that  the  fourth  integral  is  an  exact  integral,  one  obtains: 

where  X„'  =  d\ld^  . 

Similarly,  the  right-hand  side  of  Equation  9  becomes  (with  application  of  Leibnitz’  rule) 
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(10) 


[  (P  v)^  +  (pv,v)J  X^dc 


=  [pvAJ^rfC  [pv,vAJ^^/C  +  PVjvXj^^^  -  1^“  pVjVX'.JC 


C=a 


a  I-  P/i  -  9  g  ejl  rfC 

"  "  *T  I!  P.y  ■  *3  //  >•.  -  P  ^  e,  X,  ^/C 

=  -  e,  //  [p  X  rfC  -  €3  p  A  .  «3  1^“  p  a'A  -  P  «  //  ^C 

=  -  «,(,'■;  ?  1.  di)^  • f  p.,  - p  i.ij,. ., 

-  P‘i  *  '>  /.'  -PS*)  //  “'i 


Combining  Equations  10  and  11  gives 


-  //pVjVA'.rfC  +  PVjVAJ^^^ 

=  {-Pn.y  •"  ^^«(P)P.Y  -  ■"  j^-  ‘  "  ^^,(P)  +  P^,(“))«3 


(12) 


for  n=0,...,K,  where  p  and  p  are  the  pressures  on  the  top  and  bottom  surfaces  respectively,  and  P„ 
and  P„'  are  the  n*  integrated  pressures: 

P.  =  P'n  =  f^P^'.dC 


The  expression  for  the  velocity  field  in  Equation  3  is  now  inserted  into  the  left-hand  side  of 
Equation  12: 
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m»0 

K 


//p  E  KA.>-A  *  /.'p[e 


m»0 


-  //  p  E  ip.i.  E  •  E  E 

r»0 


m-O 
it  jr 


111*0  r*0 


C=a 


iii«0 

r  K 


E  P"^/,’  •  E  E  P  Kip;)  // 

— n  “  *=0  r*0  *P  • 


JT  r 


-  E  E  p  w’-'P?  /'  *  E  E  p  ip-if 

wO  fO  •  _  »  -  « 


in*0  r*0 


c=« 


E  pip„p„  *  E  E  p  K'p;),  p™  -  E  E  p  'p.'p;  yl 

«*0  ■■•O  r»0  *'  M-0  r«0 

E  E  p  'p.'p^(^.^^,)  i{:5 


r  r 


*•0  p«0 


(13) 


where 


P-  =  .  P„  =  /,’  .  y:,  =  /_’ 


(13a) 


Thus,  the  equation  for  fluid  sheets  (Equation  12)  for  an  inviscid  fluid  can  be  written  as 


E  I  p’p^mP-  *  E  p  K'P^),  p™  -  E  p  "'X p. 

*•0  I  r-O  ”  r«0 

•  E  E  p  IP.  IP?  (3.^,3,)  i':2 

M-O  r-O  “ 


(14) 


for  n=0,...,K. 
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Equation  14  can  be  reduced  further.  Prior  to  this  reduction,  some  intermediate  results 
should  be  recorded  for  later  use.  in  the  continuity  equation.  Equation  S,  the  dummy  index  m  is 
changed  to  r,  and  after  multiplying  this  equation  with  X,.  and  it  should  be  summed  over  m 
(m=0,..,K). 

This  equation  is  integrated  through  the  vertical  direction  to  yield 

I!  E  E  '/.’EE  = «  < 

— -A  — ^  _ A  _  JV 


*»0  r*0 


n«0  M) 


After  interchange  of  the  order  of  summation  and  integration  and  using  Equation  1 3a,  we  have 

E  E  Ky-.-!:  E  =  o 

r»0  ««0 


The  left-hand  side  of  Equation  1 4  is  now  considered  and  the  chain  rule  is  used  for 
differentiation  to  expand  the  second  term: 


LHS<„,  = 


E  *  E  p  -  E  p 


E  p  y-  *  E  E  p  {K^rK)  i;: 

r-O  J  ««0  r-O 


After  use  of  Equation  17  to  replace  the  third  term  in  Equation  18,  the  LHS  of  Equation  14  becomes 


=  E  p^«y^  •  E  0  x'ly^-Z  p  "'."f 

m-O  r-0  r-0 

-  E  P >'4  *  E  E  p  Ifia 

r-0  m-0  r-0  ^ 

-  E  {p»'«>’».  *  E  p  ip;  -  E  p 

*■•0  ^  r»0  r»0 

-  E  E  p 


»i«0  r-0 


The  value  of  [y„"  y„,"l  can  be  determined  by  using  integration  by  parts: 
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(19) 


[y:,  *  yj^]  =  //  -  /.’  x.a.i,)'<iC 

'  (^-^'■^.)  1(1!  -  >■"  • 


Inserting  this  identity  into  Equation  18  and  using  this  to  replace  the  left-hand  side  of  Equation  14 
yields 


jr 

E 


f*0 


pw^y-  •  i:  p  'KT  y_  *  E  p  w.ipfy, 

r»0 


(20) 


for  n=0,  ...,K. 

The  derivation  of  the  general  Green-Naghdi  equations  for  an  inviscid,  incompressible  fluid 
have  now  been  completed  .  The  governing  equations  for  inviscid  flow  are  then:  the  two  kinematic 
boundary  conditions  in  Equation  4,  the  K  conservation  of  mass  equations  in  Equation  8,  and  the  K 
approximate  conservation  of  momentum  (vector)  equations  in  Equations  14  or  20.  Note  that  the 
conservation  of  momentum  equations  are  therefore  3K  scalar  equations  for  three-dimensional  flows. 
The  variables  include  3K  unknown  components  of  W„  K  integrated  pressures  P.,  and  two  conditions 
on  the  bounding  surfaces.  On  the  top  surface  either  or  0  is  unknown,  depending  on  the  problem. 
Similarly,  on  the  bottom  surface  either  a  or  p  is  unknown.  Thus,  we  have  4K-t-2  unknowns  and  the 
same  number  of  equations,  and  the  system  is  closed. 

In  some  previous  work  and  in  our  examples,  a  so-called  "restricted"  theory  is  used.  In 
these  theories  the  X  functions  which  are  used  are  often  polynomials,  in  which  the  last  term  is 
restricted  to  have  no  components  in  the  x‘  or  x^  directions.  This  situation  will  be  discussed  separately 
below. 

We  can  make  several  observations  about  the  results  so  far.  The  equations  depend  only  on 
x‘,  x^  and  t  and  do  not  have  any  explicit  dependence  on  the  variable  ^  After  the  initial  assumption  of 
the  form  of  the  velocity  distribution  in  the  z  direction  was  made,  no  terms  were  thrown  out.  The 
governing  equations,  like  the  conservation  laws  from  which  they  were  derived,  are  Galilean  invariant. 
Because  no  scale  was  introduced,  there  is  no  explicit  flow  situation  for  which  this  theory  is  most 
applicable.  The  governing  equations  derived  this  way  are,  to  be  sure,  an  approximation.  However, 
the  limits  of  this  approximation  are  implicit  and  must  be  determined  by  numerical  or  physical 
experiment.  Even  for  the  lowest  level  theory,  the  governing  equations  are  nonlinear  because  both  the 
conservation  of  momentum  laws  and  the  boundary  conditions  are.  These  are: 
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(3) 


for  n=0,...,K. 

Because  the  governing  equations  are  approximate,  they  do  not  exactly  satisfy  Kelvin’s 
theorem  and  the  flow  computed  from  these  equations  does  not  remain  irrotation^.  Recall  that 
irrotationality  is  not  an  property  of  an  inviscid  fluid,  but  rather  a  consequence  of  an  assumption  that 
the  fluid  is  initially  quiescent  and  is  acted  upon  by  conservative  forces.  Shields  and  Webster  (1989) 
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showed  that  the  n*  level  shallow  water  fluid  sheet  theory  did  satisfy  conservation  of  circulation  in  an 
average  sense  across  the  fluid  domain  and  the  flow  does  remain  approximately  irrotational  in  an  initial 
value  problem  when  the  initial  state  was  quiescent.  However,  the  treatment  of  steady  flow  (time 
invariant)  problems  does  require  some  additional  specifications  of  the  average  circulation  (or  of  the 
vorticity  distribution). 

The  appearance  of  the  momentum  equations  in  Equation  20  is  deceptively  simple  but  note 
that  there  are  two  levels  of  implied  summation  from  index  repetition.  Actual  evaluation  of  the 
equations  is  sufficiently  tedious  that  it  is  impractical  to  carry  out  any  but  the  first  one  or  two  levels 
without  the  use  of  computer  programs  to  perform  the  calculus  and  the  algebra. 

Discussion 


The  result  of  the  derivation  above  has  been  to  reduce  the  dimensionality  of  the 
three-dimensional  equations  to  a  set  of  two  dimensional  equations  in  x',  x^,  and  t.  As  such,  these 
equations  are  reminiscent  of  equations  for  a  membrane  although,  unlike  a  membrane,  this  "fluid 
sheet"  has  a  much  greater  kinematic  complexity.  For  instance,  a  membrane  has  only  one  kinematic 
variable,  the  location  of  the  membrane  for  a  given  x‘  and  x^.  The  fluid  sheet  has  vectors  W„,  one  of 
which  may  be  identified  as  the  "location"  of  the  sheet,  but  the  others  of  which  are  clearly  kinematical 
ingredients  that  have  no  counterpart  as  a  membrane. 

Using  a  development  that  is  an  analog  of  the  development  of  the  three-dimensional 
equations  (for  instance,  the  Navier-Stokes  equations).  Green  and  Naghdi  developed  a  continuum 
model  of  two-dimensional  sheets  with  kinematic  complexity.  This  development  specifies  the  general 
form  to  be  expected  with  arbitrary  complexity.  In  this  development,  the  kinematic  ingredients  are 
called  "directors"  and  the  sheet  is  a  "directed  fluid  sheet"  or  "Cosserat”  surface. 

In  their  treatment.  Green  and  Naghdi  regard  a  set  of  governing  equations,  such  as  those 
developed  in  the  previous  section  as  a  postulated  set  of  equations  motivated  by  the  three  dimensional 
equations.  The  equations  are  to  be  validated  by  comparison  with  the  general  theory  or,  if  required, 
modified  to  reflect  the  physical  principals  embodied  in  the  general  fluid  sheet  model.  It  is  fortunate 
that  an  ideal,  incompressible  fluid  has  such  a  simple  constitutive  relation  (its  internal  stresses  are  only 
pressures  and  these  do  not  depend  on  the  rate  of  strain  of  the  fluid)  that  the  variational  procedure 
above  does  yield  a  set  of  governing  equations  that  fits  the  mold  provided  by  the  GN  fluid  sheet 
model. 


One  of  the  distinct  advantages  of  fluid  sheet  theory  is  that  it  always  results  in  approximate 
governing  equations  for  unsteady,  three-dimensional  flows.  The  specialization  of  these  equations  to 
either  two  dimensions  or  to  steady  flows  presents  no  difficulty.  However,  many  of  the  specialized 
numerical  techniques  used  in  fluid  mechanics  are  developed  with  appeal  to  specialized  methods  which 
depend  on  the  flow  being  steady  or  two-dimensional  (or  both),  and  are  therefore  limited  in  their 
applications. 


Finite  Depth  Applications 

In  this  section,  a  special  case  of  the  theory  for  shallow  water  is  provided.  In  this  case 
further  reduction  of  the  equations  is  possible  due  to  the  choice  polynomial  weighting  functions.  The 
equations  for  shallow-water  are  previously  given  by  Green,  Laws  and  Naghdi  (1974)  and  Green  and 
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Naghdi  (1984,  1986  and  1987)  and  Shields  and  Webster  (1988).  In  this  chapter  these  equations  are 
re-derived  based  on  the  general  derivation  above.  The  present  work  yields  equations  that  are  identical 
to  those  given  by  Green  and  Naghdi  (1984  and  1986);  that  is,  these  equations  can  algebraically  be 
transformed  into  those  of  Green  and  Naghdi.  In  this  research,  constraints  of  the  restricted  theory  that 
existed  in  the  earlier  work  of  Shields  and  Webster  have  been  removed,  and  a  generalized  set  of 
equations  for  an  arbitrary  level  of  the  Green  and  Naghdi  theory  is  presented.  The  complete  set  of 
equations  for  the  first  two  level  theories  for  shallow-water  applications  will  next  be  presented. 

Equations  for  Shallow  Water 

With  the  choice  of  the  polynomial  weighting  functions,  \(^)  =  it  is  possible  to  reduce 
the  general  equations  further.  In  this  case  the  various  y^,  y^,  etc.  can  be  expressed  using  a  single 
function: 


With  the  use  of  Equation  21,  one  obtains 


ymn  ~ 


(22) 


The  velocity  field  is  given  by 


x  X 

^  =  E  >  ''3  =  E 


ii«0 


fi«0 


(23) 


With  the  use  of  Equation  23,  the  equations  for  the  kinematic  boundary  conditions  and  the  conditions 
of  the  conservation  of  mass  and  momentum  may  be  obtained.  The  kinematic  boundary  conditions  are 
given  by 

X  X 

Y,  E  “.T 

f»»0 


-  E  P"  P.,  <25) 

if«0 


The  continuity  equation  becomes 


X 

E 


P"  =  P, 
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(26) 


^n.y  ^  =  0  for  n=0,l,...^-l  (29) 

Finally,  the  conditions  of  the  conservation  of  momentum  become 

E  *  E  p  "7  •  E  p 

«»0  rwO  r»0 

(-^  ^  P.r  -  “.r)  ^  -  P  «  y->  -  P  ^»(P)  ^  P 

for  n=0,...,K. 

Equations  23  through  30  are  equivalent  to  those  previously  given  by  Green  and  Naghdi 
(1984)  and  Shields  and  Webster  (1988).  It  is  noted  that  Equation  28  is  the  equation  that  is  related  to 
the  so-called  "restricted  theory."  The  details  of  the  discussion  about  the  restricted  theory  will  be 
given  in  the  following  sections. 


17 


PART  III:  IRROTATIONAL  FLOW  AND  THE  RESTRICTED  THEORY 


Discussion 

When  Green  and  Naghdi  introduced  their  theory,  they  restricted  the  last  component  of  the 
director  so  that  it  remains  vertical  at  all  times.  In  the  theory  of  Green  and  Naghdi,  specific 
constitutive  equations  are  required  for  the  3-D  response  functions;  that  is,  the  terms  on  the  right-hand 
side  of  Equation  20,  which  for  a  more  general  fluid  can  be  considerably  more  complicated.  These 
constitutive  equations  represent  the  material  properties  of  the  fluid  and  its  particular  geometry.  Also 
the  inertia  coefficients  y„  and  the  relationship  of  the  velocity  fields  v  to  the  director  velocities 
need  to  be  specified.  Green  and  Naghdi  chose  the  response  functions  so  that  the  pressure  is  the  only 
component  that  determines  the  mechanical  power.  The  constraint  responses  are  found  such  that  the 
corresponding  mechanical  power  is  zero.  They  also  assign  force  vectors  to  obtain  proper  responses. 
The  above  procedure  is  central  to  Green  and  Naghdi’s  approach  and  is  what  makes  their  theory  self- 
consistent  in  its  internal  structure.  Notably  in  the  theory  of  plates  and  shells,  their  theory  shows  its 
self-consistency  in  its  ability  to  satisfy  both  dynamic  boundary  conditions  and  kinematic  boundary 
conditions.  Many  competing  approaches  used  to  form  the  three-dimensional  equations  ended  up 
having  inconsistencies,  and  more  specifically,  both  boundary  conditions  were  not  satisfied  at  the  same 
time.  If  one  starts  with  the  right  kinematic  conditions,  one  ended  up  with  the  wrong  dynamic 
conditions,  and  vice  versa.  The  approach  of  Green  and  Naghdi  can  model  a  general  fluid  by 
specifying  its  constitutive  equation  without  any  conflict. 

Shields  (1986)  set  W^''  =  0  (corresponding  to  a  restricted  director)  in  order  to  satisfy  the 
continuity  equation.  This  was  necessary  in  order  to  obtain  K-i- 1  conditions  from  the  continuity 
equation  because  of  the  transformation  that  had  been  introduced.  Still,  no  clear  meaning  of  the 
constraint  is  offered.  It  appears  to  exclude  the  solution  which  may  be  possible  otherwise. 

The  condition  with  no  constraint  is  investigated  in  this  research.  Originally  the  restricted 
theory  meant  the  first  level  of  the  direct  theory  with  a  constrained  director.  This  concept  is  now 
extended  to  the  K*  level  theory.  It  is  called  a  restricted  theory  if  the  K'^  components  of  the 
two-dimensional  velocity  components  are  constrained  (corresponds  to  constraint  of  the  K*  component 
of  the  directors).  In  this  work  it  will  be  shown  that  this  constraint  has  a  simple  meaning. 

For  a  two-dimensional  flow  (Equation  28  in  the  previous  section),  the  most  general  solution 
of  Equation  28  becomes 

W'jf  =  constant  (31) 

Consider  now  two-dimensional,  steady  periodic  waves.  In  determining  the  wave  celerity  of 
steady  periodic  waves  an  additional  assumption  is  needed  in  order  to  ensure  that  solutions  are  unique. 
Traditionally  this  is  accomplished  by  either  of  two  definitions  of  irrotational  wave  speed  introduced  by 
Stokes.  Cokelet  (1977)  defined  the  circulation  per  unit  length,  C,  by 
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(32) 


C  =  —  f*"  u  dx 
X  Jo 

where  X  is  a  wave  length  and  u  is  a  horizontal  component  of  the  velocity.  In  the  work  of  Cokelet 
(1977),  Equation  32  is  satisfied  by  the  choice  of  reference  frame  which  travels  with  the  wave  speed  c. 
Because  the  flow  is  assumed  irrotational.  Equation  32  holds  at  every  vertical  location  in  the  fluid. 
According  to  Stokes’  theorem,  the  vertical  gradient  of  the  averaged  horizontal  velocity  is  zero  if  the 
flow  is  irrotational.  The  wave  speed  defined  by  Cokelet  is  then  that  according  to  Stokes’  first 
definition.  It  is  noted  that  Stokes’  first  definition  is  based  on  the  prior  assumption  of  an  irrotational 
flow. 


Stokes’  first  definition  has  been  used  by  most  researchers.  However,  in  the  direct  theory  no 
assumption  of  an  irrotational  flow  is  made  a  priori.  In  the  direct  theory  it  is  more  natural  to  adopt 
Stokes’  second  definition  of  the  wave  speed.  However,  to  be  consistent  with  previous  work,  Stokes’ 
first  definition  is  used  here,  and  accordingly  additional  irrotational  requirements  are  needed.  They 
can  be  treated  analogously  in  the  direct  theory.  Equation  32  will  be  used  for  the  condition  of 
irrotational ity.  We  recall  that  each  weighting  function  represents  different  vertical  dependence.  Since 
there  is  no  vertical  gradient  of  the  averaged  horizontal  velocity  for  an  irrotational  flow,  we  may 
obtain  K+ 1  conditions  for  the  requirements  of  an  irrotational  flow  if  Equation  32  is  to  be  satisfied  at 
any  vertical  location.  These  are: 

+  c)  dr  =  0  (33) 


forJ^lX-JC  (34) 


where  X  is  a  wave  length  and  c  is  the  speed  of  a  moving  frame  which  is  the  same  as  the  wave  speed. 
Equation  33  is  the  definition  of  the  wave  celerity  and  is  analogous  to  the  first  definition  of  Stokes. 
Equation  34  is  an  expression  for  global  irrotational  requirements,  and  may  be  used  for  the  measure  of 
the  vorticity.  If  Equation  34  is  not  met,  then  the  solution  is  not  an  irrotational  solution  in  these 
integral  senses.  If  Equation  34  is  satisfied,  then  Equation  33,  the  definition  of  the  wave  speed,  is 
independent  of  vertical  location  within  the  fluid  field  and  is  therefore  equivalent  to  the  first  definition 
of  Stokes. 

We  now  go  back  to  Equation  31.  In  order  to  satisfy  Equation  34,  it  is  necessary  that  W^'  = 
0.  Consider  for  example,  two-dimensional  solitary  waves.  Since  W,^'  vanishes  far  upstream  and 
downstream,  W^'  becomes  identically  zero  in  this  case.  Recall  that  W^'  =  0  corresponds  to  the 
statement  of  the  restricted  theory.  From  this  it  becomes  clear  that  the  restriction  is  an  implicit 
assumption  of  an  irrotational  flow.  This  is  considered  as  a  necessary  condition  for  the  irrotational 
flow  although  it  is  not  a  sufficient  condition.  Shields  stated  that  Theory  I  does  not  admit  shear  flow 
solutions.  This  is  not  true  in  general,  but  it  is  true  if  one  uses  the  restricted  theory.  It  is  possible  to 
model  waves  with  shear  flow  (or  current)  with  Theory  I  which  is  unrestricted.  Higher  level  theories 
are  capable  of  modeling  shear  flow  solution  with  the  restricted  theory. 
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It  becomes  clear  that  the  restricted  theory  is  needed  in  order  to  model  shallow-water  problems 
whose  fluid  field  is  considered  as  irrotational.  In  modeling  rotational  flow,  either  a  higher  level 
theory  or  the  unrestricted  first  level  theory  may  be  used  depending  on  the  accuracy  of  the  solution 
desir^  and  other  conditions,  if  necessary.  From  now  on,  the  restricted  theory  will  be  used  in 
shallow-water  problems  unless  otherwise  stated. 


Two-Dimensional  Unsteady  Flow  -  Theory  I  (General  Formulation) 

In  this  section  two-dimensional  equations  for  a  free  surface  flow  over  an  even  bottom  are 
provided.  This  section  presents  a  basic  overview  of  the  treatment  of  these  equations.  The  notation 
will  be  changed  from  tensor  notation  to  a  component  notation  in  which  x  =  x'  is  the  flow  direction. 
But  usage  of  f  is  retained  to  represent  the  vertical  direction  together  with  z  for  convenience.  The 
coordinate  system  is  taken  so  that  the  bottom  is  expressed  by  o((x,t)  =  0  with  p,  the  unknown 
pressure  on  the  bottom.  The  upper  surface  /3(x,t)  is  a  free-surface.  On  this  surface,  p  needs  to  be 
specified.  With  K=  1,  Theory  I  equations  for  this  specific  case  will  next  be  presented. 

In  Theory  I  the  velocity  profile  is  given  by 

B  =  Kq  +  M,  C  (35) 


W  =  Wj,  +  Wj  C  (36) 

where  u  is  the  horizontal  velocity  component  and  w  is  the  vertical  velocity  component.  The 
kinematic  boundary  conditions  are 

Wj,  =  0  (37) 


^1  P  =  P,  +  Ho  P, 

The  continuity  equation  yields  the  following  conditions: 

B,  =  0  (39) 

Bj^  +  W,  =  0  (40) 

Equation  39  means  that  the  restricted  theory  is  used.  The  conditions  of  the  conservation  of 
momentum  yield  the  following  equations: 
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(41) 


P  “o,  +  P  “o  “to 


2 


Ox 


(42) 


2 


1 

2 


P^  “o 


(43) 


3 


P*  «o  »*'u 


J, 

3 


(44) 


Equations  41  and  42  are  statements  of  conservation  of  the  horizontal  component  of  the  momentum 
equations  for  n=0  and  1,  respectively,  and  Equations  43  and  44  are  similar  statements  of  the  vertical 
component  of  the  momentum  equation. 


Because  there  are  many  unknown  variables  in  this  theory,  it  is  convenient  for  computational 
purposes  to  reduce  the  system  of  equations  to  one  with  fewer  unknowns  and  equations.  The  reduction 
is  usually  done  by  expressing  all  of  the  vertical  components  of  the  velocity  in  terms  of  the  horizontal 
components  of  the  velocity  through  the  continuity  equation  and  the  kinematic  boundary  conditions.  It 
is  noted  thar  this  can  be  done  for  any  level  of  the  theory.  Moreover,  some  terms  are  decoupled  from 
other  variables.  For  instance  P,  (P^  in  the  K*  level  theory)  is  decoupled  and  occurs  in  only  one 
place.  Equation  42.  Since  there  is  no  particular  interest  in  this  unknown,  it  may  be  possible  to 
discard  this  variable.  The  parameter  p  occurs  only  in  Equation  43,  and  therefore,  this  equation  is 
used  to  express  p  and  can  be  removed.  If  one  wishes,  it  is  always  possible  to  compute  p  from  the 
solutions  obtain^.  Also  Po  (Po,  Pi,...,Pk.i  for  the  K*  level  theory)  may  be  eliminated.  As  a  result, 
the  original  set  of  eight  equations  may  be  reduced  to  a  system  with  two  equations  subject  to  the 
variables  /3  and  Uq.  For  the  second  level  theory  it  is  possible  to  reduce  the  system  to  three  equations 
with  three  unknowns.  For  steady  flow  problems  one  of  the  components  of  the  horizontal  velocity 
may  be  further  reduced. 

The  number  of  unknowns  is  now  reduced.  As  was  mentioned  above,  the  vertical  components 
of  velocity  can  be  expressed  in  terms  of  the  horizontal  components  of  the  velocity  to  give 

Wo  =  0  ,  w,  =  -  «o,  (45) 


From  Equations  43  and  45,  an  expression  for  pressure  on  the  bottom  surface  is  obtained  as 
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(46) 


P 

P 


I  j  «0  “Qx* 


+  «  p 


It  is  noted  that  P,  occurs  only  in  Equation  42.  Since  there  is  no  interest  in  P,,  Equation  42  will  be 
omitted  here.  A  reduced  set  of  differential  equations  with  three  unknowns  may  thus  be  obtained. 
These  three  equations  are: 

P,  =  -  K  P], 


P  +  P  “o  “o  I  ~  P 


Px 


'Ox 

P  P 


(48) 


-  i  i  P’  •  I  p’ -  p  ^  ^ « p’ 

j  J  J  p  p  2 

Pq  may  be  eliminated  using  Equations  48  and  49  to  yield  tlie  following  two  equations  with  two 
unknowns,  /S  and  Ug: 

P.  =  -  K  P], 


ar 


dx  dxdt 


dx^dt 


dx 


dx 


iP 

dx 


.dx, 


(50) 


dUr. 


.  3  e  „  M  :u!2  -  p2r:«  u!! .  .  rj:* 

“  ftr  a,!  Sx  ”  Bx' 


Two-Dimensional  Steady  Flow  -  Theory  I  (General  Formulation) 

In  this  section,  equations  for  steady  flow  will  be  obtained  by  setting  every  time  derivative 
equal  to  zero.  Since  there  are  three  governing  equations.  Equations  47  through  49,  and  three 
unknowns  03,  Uq,  and  P,),  these  three  equations  will  be  used  here.  The  pressure  on  the  upper  surface 
is  assumed 

P  =  Pa  ■  9  (51) 


where 
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(52) 


The  variable  T  represents  the  constant  surface  tension  on  the  surface  and  p,  is  the  atmospheric 
pressure.  Since  the  fluid  is  assumed  incompressible,  p,  =  0  without  loss  of  generality  for 
convenience. 

Among  these  three  equations,  two  can  be  integrated.  Equation  50  can  be  integrated  with 
respect  to  x  to  yield 

«o  P  =  <?  (53) 

where  Q  is  a  constant  of  integration.  A  meaning  of  this  constant  can  be  easily  interpreted. 

The  mass  flux  per  unit  span  is  expressed  as 

!o  “  "  fo  “o  =  “o  P  =  <? 


Therefore  Q  represents  the  mass  flux  per  unit  span. 


Integration  of  Equations  48  and  49  yields 

P  «o“ox  =  -  — 
P 


TK 

P  ^(1^/)' 
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s*  =  \p/p  +  U^]  di  =  PJp  *  Q  Uo 


The  relation  between  S  and  S*  is  now  obvious;  S  may  be  interpreted  as  the  momentum  flux  per  unit 
span  less  the  momentum  defect  due  to  the  effect  of  surface  tension. 


With  the  use  of  Equations  54  and  58,  a  single  governing  equation  inclusive  of  /S  may  be 
obtained  from  Equation  56; 


- 1  o’  ^  ♦  2!  *  1 .  B’ _ l- 


P  P  2 


p  /(I -Pi)’ 


+  5  =  0  (60) 


Equation  60  may  be  integrated  once  more  after  multiplying  by  After  integrating,  multiply  by 

/3^  to  obtain 

I  O’  U,’  *  «  P’  •  2  «  P’  .  2  S  p  .  -  0’  ■  0  (61) 

’  Pi/i^ 

where  R  is  another  constant  of  integration.  Solutions  to  these  equations  for  different  values  of  Q,  R, 
and  S  correspond  to  regular  waves  and  solitary  waves.  Such  waves  have  been  computed  by  various 
researchers  using  these  equations  (for  instance,  see  Ertekin  1984). 


Two-Dimensional  Unsteady  Flow  -  Theory  II  (Flat  Bottom 
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Continuity  equation: 


(66) 


1/2  =  0 


(restricted  theory) 


+  wj  =  0 


+  2  Wj  =  0 


(67) 

(68) 


Conservation  of  momentum: 
n=0,  x-component : 


P  “or  I  “i,  *  I  P’  “1  “lx  I  P^  “0  “lx  +  I  P^  “1  “ox  P  “0  “( 

+  1  p2  “,  W,  +  1  p’  tt,  H-2  =  Px  ^  ^ 

2  3  P  P 


Or 


(69) 


n=l,  x-component : 

I  P'  “0.  *  I  P’  “1,  *  7  p"  “1  “u  +  I  P’  “1  “a.  »  -J  P“  “0  “o.  ♦  ^  P‘  “1 


.  1  P>  .  p  p  1  -  -i 

3  P  P 


(70) 


n=2,  x-component : 

I  P’  ^  P‘  “u  ♦  -J  P’  “1  “1,  *  7  P*  “1  “cu  *  I  P’  “0  “q.  ♦  I  P’  “1  “■j 

*  I  p* .  p-  p,  ^ 

4  P  P 


(71) 
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n=0,  f-component  : 


1 

2 


+  I  P'  W2,  +  I  P"  «,  H-i,  *  j  P’  «o  '^2*  I  P’  «i 

+  A  W,  Wj  +  -i  P*  w,*  =  -  P  +  -^  “  ^  P 


I  «0  »*’l* 


(72) 


n=l,  f-component ; 


-  P'  H-,.  *  - 
3  “  4 

2 

+  — 
5 


P"  Wj,  +  j  P*  ^  P^  \  P"  “l  "  I  P’  “O 

P  1 

p5  ^  2  +  1  p^  W,  H-,  ^P^  W,2  =  -p  p  +  -  i  ^  P^ 

^4  *  3  p2 


(73) 


n=2,  f-component : 


_1 

4 


P^  H-2,  -g  P‘  «,  M'z, 
P*  Wj*  +  I  p*  W,  ^2 


*  j  P’  «0  >^2,  *  I  P'  “l  »»'U  ^  I  P'  “O 

+  lp<  w2  =  _p2  £.  4.  2^^  -  4  «  P’ 

4  '  P  P  3 


(74) 


There  are  a  total  of  1 1  equations  to  solve.  Equations  64  to  74.  The  system  of  equations  is 
reduced  as  follows.  The  variables  Wo,  w„  and  are  eliminated  using  Equations  66-68.  The 
expressions  of  the  vertical  components  of  the  velocity  are  given  by 

1  ( 
=  0  =  “  “o »  ^2  =  “  ^  “l  X  ^ 

With  Equation  75,  five  governing  equations  for  /5,  Uo,  u,,  Po,  and  P,  are  written  here: 


P.= 


(76) 
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(77) 


P  “or  I  P*  «u  +  -g  P’  «1  “i ,  *  I  “o  “i  X  ^  p  “o  “ox 

=  ph.  .  ^ 

P  P 


|p^«0. 


I  P’  “l.  P*  “l  “l  X  ^  I  p’  “o  “J  X  *  I  P^  “o  “ox 

=  P  P,P/P  -  Pi  X  Ip 


I  P’  “Oxr  -  I  P"  “l  xr  -  ^  P*  “l  “l  XX  -  I  P^  “O  “l  «  ^  P’  “fx 

*  I  P‘  -  7  p'  - 1 P’  - 1 P’  ‘L 

-  -  p  -  »  —  - 1 « p" 

p  p  2 


p’  “ox  “l  X  -  T  P’  “l  “oxx  -  T  P^  “o  “oxx  7  P^  “o% 

10  5  4  4  (80) 

.-P=£*2^-4«p> 

P  P  3 

It  is  noted  that  Pq  and  P,  may  be  eliminated.  The  result  is  the  following  set  of  three  equations 
in  Uo,  u,.  and  |8,  and  assuming  in  this  case  that  the  surface  tension  is  zero  and  thus  ^  is  zero: 


p,  ■  -  («.  p  *  p*)^ 
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-  120^  -  «0  ,  120pi£f^2  »  60  p>^^  *  40  .  15  p>-^ 


etc  dxdt 


etc  dxdt 


120  +  120  «o-^  +  120  ^  +  20  P  (3  +  p  «,)^ 

ax  dx  dx  ax  dx 


*  ,80  P^^^^  .  60  p>#:^  -  120  PK  *  *  40 p=^^ 

dx  dx  dx  dx  dx  dx  dx^  dx  dx^  (82) 


^  -  “  "’f""  *  P  “.)  f  ^  •  30  .  12  P'|l^ 

3x  ax*  &  ax*  dx  dx*  dx  ax* 

-  10  p*(4  «o  ^  3  P  «,)^  -  3  p*(5  «o  +  4  P  11,)-^  .  and 

dx*  dx* 


-  60—  -  40  p-^  +  60p-^-?^^2  +  30  p^iP +  15 
=u  dx  dxdt 


dx  dxdt 


.6P’-^ 

dx*dr  dx*dr 


'  '““o^  *60  0^^  *5P(8«..3P.,)^ 


» •”  >’f  5  -««-•>■■>?  9  •  ■■••$9 

6  .  V,  «i...  .  »  >  SP  3'«.  a».  S'.,  ^  s..  S'., 


*  *  P  -37^  -  30  P  <“o  ♦  p  “i)  *  ‘3  P 


dx  ax* 


dx  ax* 


dx  ax* 


5  P^^ 


dx  ax* 


-  3  P*(5  «o  +  4  P  Hi)^  -  p’(6  «o  +  5  p  «j)-^ 

dx^  dx* 


Two-Dimensional  Steady  Flow  -  Theory  II  (Flat  Bottom) 

In  this  section,  the  second  level  of  the  theory  for  two-dimensional  steady  flow  with  surface 
tension  is  provided.  For  this  case,  surface  tension  is  allowed  and  the  pressure  on  the  upper  surface  is 
described  in  Equations  51  and  52.  All  of  the  time  derivatives  in  the  governing  equations.  Equations 
76  -80  are  taken  to  be  zero. 

Equation  76  may  be  integrated  to  yield 


«o  P  *  2  «i  ^ 
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As  before,  Q  may  be  interpreted  as  the  mass  flux  per  unit  span.  One  additional  term  is 
included  in  Equation  84  compared  with  Equation  53.  This  shows  the  difference  between  the  second 
level  theory  and  the  first  level  theory.  The  x-component  of  the  momentum  equation.  Equation  77  can 
be  integrated  with  the  use  of  Equation  84  to  yield 

J’o/P  •  -iVP’  -  O'l  *  -1^  *  ® 

Again,  S  may  be  interpreted  as  the  momentum  flux  per  unit  span  less  the  momentum  defect  due  to  the 
effect  of  surface  tension. 

From  Equation  72,  p  may  be  determined  if  other  variables  are  known.  As  there  is  no  interest 
in  the  variable  Pj  in  this  research,  the  corresponding  equation  is  omitted.  Then  it  is  possible  to 
eliminate  Uq  and  Pq  using  Equations  81  and  82  to  obtain  three  equations  in  three  unknowns,  /3,  u,,  and 
P,.  These  are: 


^  I,>.L  •  i  ^  P’  *  15  P’  P>.'  *  i  ^  P 

.  1  (Ji  PZ  »  — L - L  P>  pi  *  (?i  -  S  .  0 

P  P  p/ITp?  '  2  f 


^  P- ^  |<?P> ..  ♦  i<?>  e  -  -^[  P»  •  5^P‘ P‘) ^ 
*  ip’  p. «, •  j<?  P’  P.  «u  •  ip‘  P.’  *  io  P’  p/ 
-  P/  -2^  *  =  0 

4  P  3 


Further,  P,  may  be  eliminated  between  Equations  87  and  88,  but  for  these  steady  flow  equations  there 
is  no  particular  advantage  in  doing  so.  Shields  and  Webster  (1988)  present  solutions  to  these 
equations  for  both  solitary  waves  and  for  large-amplitude  regular  waves;  their  results  were  very 
satisfactory. 
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The  governing  equations  for  the  situation  of  an  unsteady  flow  over  an  uneven  bed  are 
extremely  difficult  to  derive  or  present.  As  above,  the  equations  can  be  reduced  to  three  equations  in 
Uo,  u,,  and  0.  The  governing  equations  were  obtained  using  an  algebraic  manipulation  program. 
Mathematical.  Since  these  equations  will  form  the  foundation  for  the  research  performed  herein,  a 
presentation  of  them  in  usual  mathematical  form  is  made  below.  These  equations  have  a  new 
(time-independent)  variable,  a,  the  vertical  coordinate  of  the  bottom.  The  solution  involves 
derivatives  of  a  and  it  is  assumed  that  these  derivatives  through  the  third  derivative  in  x  are  bounded. 
The  set  of  applicable  equations  so  obtained  for  an  uneven  seabed  are  listed  next  for  completeness. 
These  equations,  excluding  the  surface  tension  effects,  establish  the  basis  of  the  numerical  code  in  the 
implementation  of  mathematical  theory. 

Uneven  bottom  governing  Equation  1: 


2^- 

dt 


^  *  P  “i) 

ax  ox  ox 


-  (a  -  P)  (a  +  P)-^ 

dx 


Uneven  bottom  governing  Equation  2: 


dUn  dUn  da  dB 

-  120  (a  -  p)^  4  6-  (a  -  P)2  ^  -  120  ^  (a  -  P)-^ 

at  at  ax  ax  at 

+  120—  (a  -  P)^  — ^  -  120  (o  -  P)*  — -  +  40  (a  -  P)’  - - 

dx  ^  adx  dx  dtdx  dtdx^ 


da^  du,  du. 

^  60  ^  (a  -  p)2^  +  60  o-^  (a  -  p)^  ^  -  60  (o  -  B)(o  +  P)^ 

dx  dt  ax^  ^  ^ 

-  120  a—  (o  -  p)  +  120  o  — (a  -  p)* 

dx  dx  dt  dx  adx 


.  36  1  ^*<1 

-  60  (a  -  P)2  (a  ^  P):^  5  (a  -  p)3  (5  a  H.  3  P)  ^ 

ax  otax  dtax^ 
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120  (a  -  -  120  (a  -  p)*  +  40  («  -  p)^-^ 

cbt  ax  ax  ax  dx‘‘ 

dU  ri 

■H  120  (a  -  p)  ^  («o  ♦  ««,)  +  120  (a  -  P)  Ui  (u^  +  o  «i) 
dx  ax 


-  120 -^(o  -  p)^-^(-Uo  +  «,(«-  2p))  -  120-^(a-p)2  +  P  «i) 

dx  ax  ax  "  *  ax  ax^ 

+  120(a  -  P)2^-5^(«o  ^  P«,)  -  120-^(-  a  +  P)-^h,(«o  *  P  “,) 

dx  dx^  ^  ^ 

*  120  ^(a  -  P)-^(«o  +  ««,)  («o  ^  P«i) 
dx^  ax 


+  60  — —  JI,  +  20-^ («0  +  O  «,)l(o  -  P)*(3«o  +  (a  +  2P)«,) 
dx^  dx  ax  ax*  ax* 


,  a^u«  ,  36 

-10(a  -  P)  ^  («  ^  3p)a,)  -  60(a  -  p)*(a  +  ^p)-^-^-^ 


+  5(a  -  P)*(5o  +  3P)-^-^  +  120-^(o  -  P)-^((2a  -  P)Mo  +  «P«i)-^ 

ax*  9x  ax  ax  ax 

-  60—  +  60o-^](a  -  P)*(3«o  +  (a  +  2P)«,)-^ 

ax  ax*  j  ^ 

,  .  du. 

*  20(o  -  P)  (3(0  +  P)  «o  +  (4  o*  +  o  P  +  P*)«,)  — 

OX 
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-  60(a  -  p)2(«  +  +  10(0  -  P)’(o  +  3P)-^-^ 

dx  ax  ax  dx^ 

^  60(0  -  p)2(o  +  P)-^(«o  ^  P«,)-^ 

dx  dx^ 


-  40  -^(o  -  P)2((4o  -  p)0o  *  «(«  ^  2p)«,)-^ 

dx  3x^ 

,  -  ,  d«,  d^tt, 

+  (o  -  p)^(7o*  +  21oP  +  12p*)--i - i 

&  dx^ 

-  (o  -  P)^5(5o  >  3p)«o  M7  o^  +  21  oP  +  12  P*)«,)-^ 

dx* 


Uneven  Bonom  Governing  Equation  3: 

-  60  (o  -  P)(o  ^  P)^  ^  20  (o  -  P)*  (2o  ^  P)^ 

OT  dx*  c» 

(91) 

-  60-^  (o  -  P)(o  +  P)  >  40  (o  -  P)2  (2o  +  P)-^ 

dx  dx  dt  dx  dtdx 

-  60  (o  -  P)*  (o  *  P)-^  ^  +  5(0  -  P)’  (5o  +  3P)-^ 

dx  otax  dfdx* 

+  20-^  (o  -  P)*(2o  +  p)  ^  ^  20o  (o  -  p)*  (2a  +  P)-^ 
ax  ot  dx^  at 


32 


-  40  (o  -  p)(a^  +  aP  +  P^)-^  -  60a  (a  -  p)  (a  +  P)-^-^ 

dt  dx  ax  at 

+  40o-^  (a  -  p)"(2a  +  p)  ^  -  30  (a  -  P)^  (a  +  p)^^^ 
dx  dtdx  dx  dtdx 

+  2  (a  -  p)^  (8o*  +  9oP  +  3p*) - = 

dtdx^ 


60(a  -  P)(a  +  P)-^g  -  60(a  -  p)^a  P)4r-^ 
dx  dx  dx 

,  dUa  3^«a  3Ka 

+  5(0  -  p)’  (5o  +  3P)^^  +  60  (o  -  p)(a  +  P)  ^  («o  ^  «“i) 

dx  dx^  ^ 


+  60— (O  -  P)(o  +  P)  U^  (Uq  +  OMj) 
ax 


(.g  r'a  ^0  1C  da  df^a  ,  d^a  .  v 

'I** - -  +  15 - u,  +  5 -  («o  +  au.)  . 

\  dx^  dx  dx  dx^  ^  dx^  °  M 


(o  -  p)*(4(2a  +  p)«o  +  3(0  +  P)*«,) 


_  60  (o  -  P)  (o  ^  P)  -^  -^  (-  u,  -  2P)«,) 

dx  dxdx 

^  60(0  -  P)2(o  +  P)^  ^  («o  -  P«,)  +  60  (o  +  P)  (o  -  P)-^«,(«o  P«i) 

dx  dx^  dx  ox 

+  60-^(o  -  p)(o  +  P)-^(«o  *  ““i)(“o  *  P“i) 
dx^  dx 

-  5-^(0  -  (8(20  +  P)Uo  +  (o^  +  14op  +  9p")«,) 

dx  dx^ 

-  (a  -  P)’  ^  (5(5o  ^  3p)Uo  ^  (7a^  +  21oP  +  12p2)«,) 
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-  30(0  -  (0  +  P)  (o  +  3p)-^-^-^  +  2(0  -  P)^  (8o*  +  9oP  +  3p")-^-^ 

dx  ax  ax  dx^  or 

-  15^  (a  -  p)2  (4  (2o  +  P)  «o  +  3(0  +  P)"  «») -^ 

dx  ax 

-  ISo-^  (o  -  P)*  (4  (2o  +  P)  uo  +  3(0  +  p)"  K,)-^ 

cbt^  dr 

*  60^  (o  -  p)  (o  +  P)  -^((20  -  P)«o  ^  oP«,)-^ 


+  5  (o-P)  (8  (o^+op  +  p^)  tt(,+3  (o  +  p)(3  o^+P*)  “i)-^ 

-  30(«  -  P)'(a  >  P)'  P-^^  +  («  -  P)’(7o'  +  21aP  +  12P')^^ 

dx  dx  ax  dx^ 

d^« 

-  15— (o  -  p)^  (ao^  +  2oP  -  p^)  «o  +  20(0  +  p)^«i) 

dr  dr^ 


^  30(0  -  P)^(o  +  p)2  |&(«o  +  p«,)-^7  +  5(0  -  P)^(o  +  P)^ 

dr  dx^ 

d^u. 

3  „  ■»  _ L 


-  (o  -  p)’  (2(80*  +  9oP  +  3P0  Uo  +  5(0  +  p^u,) 


dr^ 


d«j  d*«, 


34 


PART  IV:  SOLUTION  ALGORITHM 


Integration  of  evolution  equations,  such  as  Equations  81  -  83  is  now  considered.  These 
equations  are  very  complex  from  an  algebraic  standpoint,  but  in  fact  may  be  integrated  with  little 
difficulty.  Together  they  form  a  system  of  three  coupled,  partial  differential  equations  that  are 
first-order  in  time  and  third-order  in  space.  Further,  they  generally  have  boundary  conditions  at  both 
ends  of  the  domain  (two-point  boundary  conditions).  The  key  to  an  efficient  scheme  for  their  solution 
lies  in  the  fact  that  the  highest-order  mixed  derivatives  are  only  first  order  in  time  and  second-order 
in  space.  The  scheme  presented  below  has  been  developed  to  handle  similar  problems  and  has  been 
used  extensively  to  integrate  evolution  equations  for  deep-water  waves  (Webster  and  Kim  1990). 

There  is  no  specific  limit  to  three  governing  equations;  in  fact,  any  number  is  allowable,  so  long  as 
there  are  sufficient  boundary  conditions. 

Ertekin  (1984)  devised  a  scheme  for  solving  similar  equations  and  boundary  conditions 
corresponding  to  Theory  I.  His  analysis  was  simplified  since  there  were  only  two  governing 
equations,  and  these  were  not  coupled  in  the  time  derivatives  of  the  variables.  The  equations  are 
identical  to  Equation  50.  The  first  equation  could  be  integrated  directly.  The  second  equation, 
implicit  in  the  spatial  coordinate,  was  expressed  as  a  tri-diagonal  system  of  linear  equations  and 
solved  using  the  Thomas  algorithm. 

The  Theory  II  equations.  Equations  81  -  (83,  are  considerably  more  complex  than  those  of 
Theory  I.  Moreover,  the  last  two  equations.  Equations  82  and  83,  have  the  disadvantage  of  being 
coupled  in  the  time  derivatives  of  the  variables.  However,  it  is  still  possible  to  use  essentially  the 
same  efficient  scheme  for  their  solution.  The  current  procedure  results  from  a  hybridization  of 
Ertekin’s  method  and  the  scheme  devised  by  Newman  (1968)  for  coupled  ordinary  differential 
equations. 

In  the  interest  of  future  work  with  the  theory,  this  algorithm  will  herein  be  described  in 
general  terms.  In  particular,  consideration  will  not  be  limited  to  the  three  equations  at  hand,  but 
instead  a  more  general  system  of  equations  will  be  considered  which  would  be  applicable  to  a  higher 
order  approximation  of  the  theory. 

A  system  of  K  coupled,  quasi-linear  partial  differential  equations  will  be  considered  in  the  K- 
dependent  variables.  The  variables  are  expressed  as  a  K-dimensioned  vector,  ^(x,  t),  and  the 
equations  have  the  special  form: 


A  i  ^  B  *  C  =  g  (92) 

where  A,  B,  and  C  are  K  x  K  matrices,  g  is  a  K-dimensioned  vector.  A,  B,  C,  and  g  are  perhaps 
functions  of  x  and  |  and  its  spatial  derivatives,  although  this  dependence  will  not  be  shown  in  the 
interest  of  simplicity.  The  dot  over  ^  signifies  a  derivative  with  respect  to  time.  It  is  assumed  that 
the  problem  is  posed  as  a  two-point  boundary-value  problem  in  x  and  an  initial-value  problem  in  t. 

The  domain  of  x  over  which  a  solution  to  the  equations  is  desired  is  assumed  to  be  a  uniform 
grid  of  x’s  spaced  a  distance  Ax  apart.  The  i-th  point  on  the  grid  will  be  denoted  by  x^  =  i  Ax,  i  = 
l,ns.  Time  is  also  assumed  to  be  discretized  with  intervals  At,  with  tj  =  j  At.  The  value  of  the 
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solution  vector  f  (Xj,^)  will  be  denoted  by  and  similar  superscripts  will  be  used  for  the  other 
vectors  and  matrices.  The  spatial  derivatives  will  be  approximated  by  central  differences 


(93) 


With  these  approximations.  Equation  92  can  be  written  as 

^  gtA  ^  jjtw  ^  giiA  (94) 


where 


a7  "  2Ax 


-2  —  + 
Ax* 


(95) 


Ax*  *  2Ax  J 


Suppose  that  the  solution  ^  is  known  as  the  result  of  a  boundary  condition  at  this 

point.  Then  the  solution  at  i  =  1  can  be  readily  found  from  Equation  93  as 


V*  =  -  a"'’  d  - 


(96) 


where 


(97) 


This  process  is  continued  throughout  the  domain.  That  is. 
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(98) 


At  the  boundary  i  =  ns  it  is  assumed  that  is  known  and  given  as  a  boundary  condition  at  this 

point.  Then 


e(«j) 


(99) 


With  the  value  of  thus  determined,  the  other  values  of  can  be  determined  by  back 
substitution 


^  ^(>V)  (100) 

The  values  of  are  then  used  to  estimate  the  values  of  (i.e.  at  the  next  time  step)  by 

^(iV*0  =  ^OV)  +  (101) 

This  estimate  is  only  first-order  accurate  and  would  be  unsatisfactory.  However,  we  can  obtain  an 
estimate  for  by  using  these  new  values  and  reapplying  the  procedure  at  t  =  (j-H)  At.  A  new 
estimate  for  can  be  formed  by 

^  (102) 


This  new  estimate  is  now  second-order  accurate  in  both  At  and  Ax. 
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PART  V:  NUMERICAL  IMPLEMENTATION 


A  separate  report  in  this  series  will  present  a  detailed  description  of  the  Fortran  program 
based  on  the  shallow-water  Level  II  Green-Naghdi  theory  developed  herein.  This  program  was 
developed  primarily  for  the  use  of  Corps  of  Engineers  projects  associated  with  military  and  Civil 
Works  and  is  far  superior  in  performance  and  in  programming  style  to  the  programs  used  by  Shields 
(1986)  and  Shields  and  Webster  (1988)  for  the  previous  shallow-water  work.  The  program  was 
specifically  written  to  be  generally  applicable  to  two-dimensional  solutions  of  all  sorts  of  GN  theory 
water  wave  problems. 

The  flow  chart  included  herein  (Figure  1)  describes  the  general  modular  structure  of  this 
numerical  model.  The  part  of  the  program  which  is  unique  to  the  particular  level  of  the  theory  (and 
type  of  GN  theory  for  deep-  or  sh^low-water)  is  the  subroutine  called  coeff(J,neq).  The  coeff 
subroutine  in  the  program  corresponds  to  Equations  89-91  of  this  report  (Level  II  shallow-water 
evolution  equations);  this  subroutine  is  complicated  and  critical.  The  program  was  checked  for 
accuracy,  and  preliminary  computations  using  it  have  been  performed. 

The  program  consists  of  a  main  routine  which  provides  for  the  input,  output  and  the  flow  of 
the  information.  The  principal  subroutines  are  solve(neq),  a  coding  of  the  Thomas  algorithm 
presented  above,  coeff(J,neq)  which  includes  the  GN  equations,  invmat(a,b,neq,n),  a  standard  linear 
equation  inversion  routine  using  Gauss-Jordan  elimination,  and  filter,  a  digital  smoothing  filter  to 
remove  the  spurious  ripples  near  the  wavemaker.  A  routine  bottom  determines  the  various  spatial 
derivatives  of  the  bottom  topography  required  in  the  theory.  Several  small  routines  to  perform 
standard  vector  and  matrix  operations  are  included  to  make  the  program  self-sufficient. 

For  the  shallow-water  wave  study,  it  is  assumed  that  the  time  history  of  the  waves  is  known 
at  the  left-hand  boundary  of  the  domain.  The  waves  are  input  not  only  as  a  local  wave  height  history, 
/3(t)  but  also  as  a  history  of  the  corresponding  values  of  the  other  variables  in  this  Level  II  theory 
[Uo(t)  and  u,(t)l.  These  variables  are  obtained  from  the  solution  to  the  steady  flow  equations 
(linearized  for  small  wave  amplitude)  and  only  roughly  represent  the  flow  which  corresponds  to  steep 
waves.  These  linear  solutions  are  (for  waves  proceeding  from  left  to  right  with  a  celerity  c): 


P(0  =  ~  "  ^01 


«o  =  Po 
“i  "  Po 


12g[20  .  nm _ 


c[240  +  mikhf  +  3(JfcA)^] 

nogjkhf 

cA[240  +  m(khf  +  dikh)*] 


cos[(i)r] 


c  = 


24gh[{khf  *  10] 

\  [240  +  104 (khf  +  3ikJi)*] 


(103) 


and  where  jSo  is  the  wave  amplitude  and  h  is  the  water  depth.  In  that  these  linear,  small-amplitude 
waves  are  not  exact  boundary  conditions  for  finite  amplitude  waves,  there  are  some  small  oscillations 
near  the  wave  maker  as  the  solution  "finds  itself."  This  is  similar  to  the  flow  near  a  flap-type  wave 


38 


maker  in  a  physical  wave  tank,  since  the  boundary  condition  at  the  flap  is  only  roughly  like  a  free 
water  wave.  These  numerical  ripples  cause  problems  if  allowed  to  propagate  and  are  thus  filtered  out 
near  the  wave  maker.  Free  waves  elsewhere  are  not  filtered.  The  historical  wavemaker  input  values 
are  placed  sequentially  in  the  first  three  elements  of  the  variable  arrays,  so  that  spatial  derivatives  can 
be  formed.  An  Orlansky  free  boundary  condition  is  placed  at  the  other  end  of  the  domain. 

The  equations  are  integrated  and  the  results  are  dumped  into  several  files  as  they  become 
available.  Appropriate  output  for  use  at  CERC  may  be  tailored  to  the  desired  user  needs  in  a  given 
application.  An  important  feature  of  this  program  is  that  it  has  re-start  capability  and  the  program  is 
arranged  to  be  able  to  continue  using  previous  computations.  This  feature  will  be  particularly  good  if 
one  is  using  a  microcomputer. 
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PART  VI:  CONCLUDING  REMARKS 


This  report  details  the  development  and  some  of  the  philosophy  behind  the  Green-Naghdi 
theory  of  fluid  sheets.  The  fundamental  principals  of  the  Level  I  and  Level  II  Green-Naghdi  theories 
are  described  for  completeness.  Derivation  of  both  the  Level  I  and  Level  II  theories  have  been 
presented  in  a  systematic  manner  due  to  complexities  of  these  theories  and  high  level  of  mathematics. 
The  powerful  symbolic  manipulator  Mathematical  was  used  for  the  formulation  and  derivation  of 
equations.  A  modular  Fortran  program  has  been  developed  for  shallow-water  Level  I  and  II  theories. 
The  2-D  numerical  model  developed  in  this  research  is  applicable  to  solutions  of  all  Green-Naghdi 
theory  water  wave  problems.  The  program  allows  a  general  description  of  the  bottom  contours  and 
of  boundary  conditions  at  the  right-hand  side  of  the  domain  (either  open  or  reflective). 

It  is  recommended  that  the  Level  III  Green-Naghdi  theory  be  developed  in  the  next  phase  of  this 
research  for  coastal  applications.  The  Level  I  and  II  theories  presented  herein  can  be  easily  extended 
to  3-D  flows,  although  the  equations  will  be  considerably  more  complex.  The  3-D  Green-Naghdi 
equations  may  for  all  practical  purposes  be  beyond  human  capabilities,  but  these  are  easily 
manageable  with  the  use  of  mathematical  and  symbolic  manipulators  such  as  Mathematica™. 

A  3-D  Green-Naghdi  theory  will  be  better  suited  for  military  applications,  including  the  Logistics- 
Over-The-Shore  operations  and  coastal  wave  problems,  since  it  can  precisely  describe  effects  of 
amplitude  nonlinearities,  frequency  dispersion,  refraction,  shoaling,  reflection,  and  diffraction  of 
waves  propagating  over  arbitrary  water  depths,  seabed  topography,  and  non-uniform  boundaries.  The 
problems  of  wave  breaking  and  re-formation  may  realistically  be  represented 
dependent  Green-Naghdi  theory.  The  solution  algorithm  for  3-D  dimensional 
redeveloped  to  accommodate  the  new  domain  and  several  design  problems  in 
both  for  military  and  civil  works. 


only  by  a  3-D  time- 
flows  will  have  to  be 
coastal  applications. 
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